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The  adsorption  of  dimethyl  methylphosphonate  (DMMP)  on  the  (0  1  0)  surface  of  anatase  Ti02,  which  is 
isostructural  with  the  (1  0  0),  has  been  studied  using  density  functional  theory  and  two-dimensionally- 
periodic  slab  models.  The  experimentally-observed  faceting  of  this  surface  has,  for  the  first  time,  been 
included  in  the  modeling.  The  relaxations  of  bare  surfaces  both  with  and  without  faceting  are  similar, 
leading  to  an  atomic-scale  roughening  due  to  inward  (outward)  displacement  of  fivefold-coordinated 
Tisc  (sixfold-coordinated  Tec)  atoms  together  with  outward  displacement  of  threefold-coordinated  O3C 
atoms.  Molecular  adsorption  occurs  by  formation  of  a  Tisc-  ■  0=?  dative  bond  with  one  or  more 
C— H-  ■  ■O2C  bonds  between  CH3  groups  and  unsaturated,  twofold-coordinated  (02c)  sites.  The  energies 
for  molecular  adsorption,  obtained  using  the  B3LYP  functional,  are  virtually  identical  (about 
-21.0  kcal/mol)  for  the  two  surfaces  and  are  also  close  to  those  found  elsewhere  for  the  rutile  (1  1  0) 
and  anatase  (101)  surfaces.  A  possible  first  step  in  the  dissociative  adsorption  of  DMMP  has  also  been 
modeled  and  is  found  to  be  thermodynamically  favored  over  molecular  adsorption  to  a  degree  which 
depends  on  faceting. 

Published  by  Elsevier  B.V. 


1.  Introduction 

The  surface  structure  and  properties  of  Ti02,  in  both  the  rutile  and 
anatase  forms,  have  been  intensively  studied  both  experimentally 
and  computationally  [1],  For  either  form,  the  adsorption  of  a  wide 
variety  of  molecules  has  been  investigated  on  several  different  crys¬ 
tallographic  planes.  An  especially  important  aspect  is  adsorption  on 
the  isostructural  (10  0)  and  (0  1  0)  surfaces  of  anatase.’  Many  of  the 
technological  applications  for  Ti02.  in  which  the  material  is  used  as  a 
catalyst,  a  photocatalyst  [2]  or  a  catalyst  support,  involve  high-sur- 
face-area  (HSA)  powders  prepared  by  high-temperature  oxidation  of 
TiCU.  The  Degussa  P25  powder,  which  has  been  widely  used  in  sur¬ 
face-chemistry  experiments  on  Ti02,  is  an  example  of  such  material. 
Typically  these  powders  are  ~80%  anatase,  and  experiment  [3-5] 
shows  that  in  this  case  the  (0  1  0)  is  the  dominant  surface.  The 
(1  0  0)  is  also  known  to  make  a  major  contribution  to  the  surface  area 
of  anatase  nanocrystals  [6]  and  to  be  especially  active  in  photocataly¬ 
sis  [7],  Recent  work  by  Dzwigaj  et  al.  [8]  has  clearly  shown  that  the 
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’  The  designations  (1  0  0)  and  (0  1  0)  will  be  used  interchangeably  since  the  two 
anatase  surfaces  are  isostructural,  as  are  the  (101)  and  (Oil).  Some  of  the  literature 
on  anatase  powder  refers  to  the  (0  1  0),  and  that  labeling  is  used  here  except,  in 
Section  1,  when  citing  works  which  refer  explicitly  to  (1  0  0). 
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morphology  of  anatase  powder  (i.e.,  the  distribution  of  exposed  crystal 
planes)  depends  on  the  method  of  fabrication  and  treatment.  How¬ 
ever,  in  general  the  (10  0)  surface  appears  to  have  a  significant,  if 
not  always  dominant,  presence  in  many  forms  of  anatase  powder. 

Experimental  data  for  this  surface  are  scarce,  due  to  the  diffi¬ 
culty  in  obtaining  anatase  samples  suitable  for  single-crystal  stud¬ 
ies,  However,  several  computational  studies  of  adsorption  have 
been  reported,  focusing  on  atomic  H  [9-11],  NH3  and  NO  [12], 
CO2  [13],  H2O  [12,14],  O2  [14],  chlorobenzene  [15],  chlorophenol 
[16],  N  heterocycles  [17],  acetic  acid  [18]  and  oxalic  acid  [19].  In 
most  cases  unsaturated,  fivefold-coordinated  (Tisc)  sites  are  active 
as  Lewis-acid  adsorption  sites.  In  particular  the  adsorption  behav¬ 
ior  of  the  bidentate  reagent  oxalic  acid  on  anatase  (10  0)  differs 
from  that  on  rutile  (110)  due  the  differing  local  arrangement  of 
Tisc  surface  sites  [19].  Unsaturated,  twofold-coordinated  (02c)  sites 
can  also  be  involved  in  hydrogen  bond  (H-bond)  formation.  The  ef¬ 
fects,  on  the  structure  of  the  (1  0  0)  and  other  anatase  surfaces,  of 
hydroxylation  and  of  steady-state  ambient  pressures  of  H2,  H2O 
and  H2S  have  been  thoroughly  investigated  in  a  series  of  density 
functional  theoiy  (DFT)  studies  by  Arrouvel  et  al,  [20,21].  These 
studies  have  provided  important  insight  into  the  properties  of  real 
Ti02  catalysts. 

Fig.  1  shows  schematic  models  for  the  (0  1  0)  surface  which,  ex¬ 
cept  for  the  (1  01),  has  the  lowest  surface  energy  (cr)  for  anatase 
[20,22],  Fig,  la  shows  the  ideally-terminated  (unrelaxed)  surface. 
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Fig.  1.  The  anatase  (0  1  0)  surface  with  the  bulk  unit  cell  shown  by  the  dashed  lines. 
Blue  (red)  spheres  are  Ti  (0).  (a)  Shows  the  ideally-terminated  surface  and  (b) 
shows  a  (1  X  2)  faceted  structure  [23],  Both  structures  are  shown  before  relaxation. 
The  different  models  are  not  to  scale.  The  encircled  section  in  (a)  shows  what  is 
removed  in  forming  the  facet.  Various  surface  sites  are  labeled  (see  text),  and  the 
region  referred  to  as  the  “channel”  in  the  text  is  indicated.  The  plane  shown  in  (b)  is 
the  (Oil),  which  is  isostructural  with  the  (10  1).  (For  interpretation  of  the 
references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of 
this  article.) 


and  Fig.  lb  shows  the  (1  x  2)  faceted  surface  (also  unrelaxed) 
obseiT/ed  [23]  in  scanning  tunneling  microscopy.  Faceting  of  the 
(10  0)  surface  exposes  the  more  stable  (1  01),  which  is  isostruc¬ 
tural  with  the  (0  1  1).  Computation  [20,22]  shows  a  =  0.44  and 
0.53  J/m^  for  the  relaxed  (10  1)  and  (10  0)  surfaces,  respectively, 
vs.  0.31  J/m^  for  rutile  (1  10).  The  effect  of  faceting  on  the  anatase 
(1  0  0)  (7  is  unknown  at  present.  All  (1  0  0)  surface  Ti  sites  are  Tisc 
while  half  the  surface  0  atoms  are  O2C  and  half  are  saturated,  three- 
fold-coordinated  O3C  sites.  With  equal  numbers  of  singly-unsatu- 
rated  Tisc  and  O2C  sites  the  surface  is  auto-compensated,  meaning 
that  the  excess  electron  density  on  a  Tisc  is  transferred  to  an  O2C 
to  maintain  a  closed-shell  configuration  on  both.  Hence,  there  is 
no  indication  [10,17]  of  surface  states  in  the  band  gap  of  the  re¬ 
laxed  but  unfaceted  surface.  Adsorption  on  the  faceted  surface, 
which  has  not  to  our  knowledge  been  studied  previously,  should 
be  especially  interesting.  For  an  appropriately-chosen  molecule 
the  structure  of  this  surface,  with  rows  of  Tisc  sites  set  below  and 
in  close  proximity  to  rows  of  O2C  sites,  holds  out  the  possibility 
of  both  dative  bonding  to  a  Tisc  and  H-bonding  to  multiple  O2C 
sites. 

The  purpose  of  this  work  is  to  examine  the  adsorption  of  a  fairly 
complex  molecule,  (CH30)2(CH3)P=0  (dimethyl  methylphospho- 
nate,  DMMP),  on  the  anatase  (0  1  0)  surface  both  with  and  without 
faceting.  Our  interest  in  DMMP  stems  from  the  use  of  this  species 
as  a  simulant  for  chemical  warfare  agents  (CWAs)  such  as  Sarin 
and  the  demonstrated  utility  of  Ti02  in  the  catalytic  destructions 
of  CWAs  (Refs.  [24-28]  and  works  cited).  The  growing  interest  in 
HSA  Ti02  powders  for  this  application,  and  the  wide-spread  use 
of  DMMP  in  experimental  work  relating  to  CWAs,  makes  such  a 
study  timely.  For  two  reasons,  the  present  work  focuses  on  adsorp¬ 
tion  of  DMMP  on  bare  surfaces,  free  of  defects  or  co-adsorbed  spe¬ 
cies.  Firstly,  the  adsorption  of  DMMP  on  Ti02,  even  under  these 


ideal  conditions,  remains  controversial.  Secondly,  the  experimental 
data  for  this  system  [29-32]  with  which  we  will  compare  our  re¬ 
sults  were  all  obtained  in  vacuum  or  in  ultra-high  vacuum  (UHV). 


2.  Computational  and  modeling  details 

2.3.  Geometry  optimization  and  energy  calculations 

All  calculations  were  done  using  two-dimensionally-periodic 
slab  (2-DPS)  models  and  the  CRYSTAL  06  code  [33-35]  which  em¬ 
ploys  Bloch  functions  constructed  from  localized  Gaussian  basis 
sets.  To  reduce  the  computational  cost,  geometry  optimization 
was  done  at  the  restricted  Hartree  Fock  (RHF)  level.  This  has  pre¬ 
viously  been  shown  [36,37]  to  give  results  for  the  bulk  and  sur¬ 
face  structures  of  Ti02  which  are  very  close  to  those  obtained 
using  DFT.  Optimization  at  the  RHF  level  has  also  been  shown 
to  give  reliable  results  for  molecular  geometries  [38[.  For  Ti,  the 
Hay-Wadt  small-core  effective-core  pseudopotential  (ECP)  was 
used  with  the  411(3nd)G  basis  set  previously  optimized  for 
perovskite  titanates  [33,39].  For  0,  a  6-31G(d)  all-electron  basis 
set  [40]  was  used  with  the  exponents  of  the  outermost  shells 
re-optimized^  for  bulk  rutile.  The  optimized  parameters  were 
0.27281  and  0.31201  Bohr^^,  respectively,  for  the  sp  and  d  shells. 
The  bulk  lattice  constants  and  atom  positions  optimized  using  this 
procedure  show  good  agreement  with  experiment  (Table  1 ).  In  sub¬ 
sequent  geometiy  optimizations,  the  lattice  constants  were  fixed  at 
these  values  while  the  positions  of  all  Ti  and  0  atoms  (as  well  as 
those  of  the  DMMP  atoms)  were  allowed  to  vary.  For  DMMP,  6- 
31G(d,p)  basis  sets  from  a  standard  database  [40]  were  used  with¬ 
out  modification. 

After  geometry  optimization,  a  single-point  calculation  of  the 
total  energy  was  done  using  DFT  with  all-electron  basis  sets  and 
the  B3LYP  functional,  which  has  been  shown  [37,41]  to  give  reli¬ 
able  results  for  the  bare  surfaces  of  both  rutile  and  anatase.  A  dis¬ 
cussion  of  the  use  of  B3LYP  in  describing  the  adsorption  of  DMMP 
is  given  elsewhere  [42  [.  For  Ti02  the  basis  sets  are  designated  86- 
411(d31)G  for  Ti  and  8-411dlG  for  O,  while  the  6-31G(d,p)  basis 
sets  were  again  used  for  DMMP.  For  Ti  and  0  the  outermost  sp 
and  d  shells  have  been  optimized  [41  [  for  Ti02,  and  the  complete 
basis  sets  are  obtained  by  combining  these  with  the  inner  shells  gi¬ 
ven  elsewhere  [33[.  Preliminary  calculations  using  somewhat  lar¬ 
ger  6-311G(d,p)  basis  sets  [40]  for  DMMP  led  to  difficulties  in 
converging  the  self-consistent  field  (SCF).  This  may  result  from 
the  small  exponent  in  the  outermost  sp-shell  of  the  P  atom 
(0.0684930  Bohr"^)  which  can  cause  linear-dependency  problems 
in  calculations  for  periodic  structures  [33[. 

The  bulk  calculations  used  an  (8  x  8  x  4)  sampling  grid,  giving  a 
k-point  spacing  of  ~0.1  atomic  units  (au).  Calculations  for  the 
(4x1)  and  (4  x  2)  supercells  used  in  adsorption  (see  Section  2.2) 
were  done  with  (1x2)  and  (1x1)  grids  respectively,  giving  a  k- 
point  spacing  of  ~0.2  au.  Reducing  the  k-point  spacing  to  0.1  au 
for  the  (4x1)  supercell  had  a  negligible  effect  on  the  energy.  Fol¬ 
lowing  Labat  et  al.  [37],  an  “extra-large”  (75,974)  integration  grid 
was  used  for  all  DFT  calculations.  In  CRYSTAL,  truncation  of  the 
sums  of  Coulomb  and  exchange  terms  in  the  Fock  matrix  is  deter¬ 
mined  by  five  overlap  criteria  [33,35]  (T1-T5).  These  were  set  at 
10“^  forTl-T4  and  10“^^  forT5.  A  high  numerical  accuracy  in  inte¬ 
gration  is  generally  required  forTi02  [37,41]  and  is  further  neces¬ 
sitated  by  the  small  exponent  (0.09983  Bohr^^)  in  the  outermost 
sp-shell  of  the  P  6-31G(d)  basis  set.  In  the  multipolar  expansion 
zone  [33-35]  a  maximum  order  of  L  =  6  was  used.  The  geometiy 
convergence  criteria  were  1  x  10“^  Hartree  for  energy,  1.8  x  10“^ 


^  Basis-set  optimization  was  done  using  CRYSTAL  03  with  the  LoptCC  script  (C. 
Zicovich-Wilson;  Universidad  Autonoma  del  Estado  de  Morelos;  Cuernavaca,  Mexico). 
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Table  1 

Observed  and  calculated  Ti02  lattice  parameters.® 


Rutile'® 

Anatase*^ 

Expt. 

Calc. 

Expt. 

Calc. 

a 

4.5937 

4.564 

3.7845 

3.781 

c 

2.9587 

2.987 

9.5143 

9.536 

u 

0.30478 

0.3060 

0.20806 

0.2063 

Volume 

62.435 

62.224 

136.268 

136.355 

®  Experimental  values  are  from  neutron  diffraction  data,  Ref.  [46].  The  lattice 
constants  (a,  c)  are  in  Angstroms,  and  u  is  a  fraction  of  a  lattice  constant.  The  vol¬ 
ume  of  the  crystallographic  unit  cell  (which  comprises  two  primitive  cells  in  the 
case  of  anatase)  is  given  in  fl?. 

The  rutile  Ti  positions  are  at  (0, 0, 0)  and  (I/2, 1/2,  Vi).  The  O  positions  are  at 
±(u,  u,  0),  (-U  +  'A  u  + 14, 14)  and  (u  + 14,  -u  + 14, 14). 

'  The  anatase  Ti  positions  are  at  ((0, 0, 0),  (14, 14, 14))  +  ((0, 0, 0),  (0, 14,  V4)).  The  0 
positions  are  at  ((0,0,0),  (14, 14, 14))  +  ((0, 0,  u),  (0, 14,  U  +  V4),  (14, 0, -u  +  V4), 
(V2, 14,  -u  + 14)). 

Bohr  for  maximum  displacement  and  4.5  x  Hartree/Bohr  for 
maximum  gradient. 

2.2.  Ti02  model  structures 

The  2-DPS  model  for  the  unfaceted  surface  used  a  (4  x  1 )  super¬ 
cell  with  six  layers,  giving  48  Ti02  units  per  cell.  The  lateral  dimen¬ 
sions  were  chosen  to  give  a  low  coverage  of  adsorbed  DMMP  (thus 
avoiding  intermolecular  interaction)  while  still  being  computa¬ 
tionally  affordable.  The  surface  area  of  the  supercell,  containing 
one  adsorbed  DMMP,  was  1 5.1 3  A  x  9.54  A  =  1 44.2  A^,  and  the  cov¬ 
erage  was  1 /8-monolayer  (one  DMMP  per  eight  surface  Tisc  sites). 
A  slab  thickness  of  six  layers  has  been  shown  [22,37]  to  be  ade¬ 
quate  for  convergence  of  the  surface  energy  of  anatase  (1  0  0). 
For  the  faceted  surface  the  supercell  was  increased  to  (4  x  2)  in  or¬ 
der  to  accommodate  the  facet.  After  removal  of  the  appropriate 
atoms  (cf.  Fig.  1 ),  the  supercell  comprised  80  Ti02  units.  In  forming 
the  facet,  equivalent  atoms  on  either  side  of  the  slab  were  removed 
so  as  to  maintain  inversion  symmetry.  Complete  relaxation  of  the 
faceted  surface  will  be  seen  (Section  3.1)  to  require  a  thicker  slab. 
Thus  calculations  for  the  bare  faceted  surface  were  also  done,  for 
comparison,  for  a  thickness  of  12-layers  and  a  (1  x  2)  supercell. 
Since  localized,  rather  than  plane-wave,  basis  sets  are  used  one 
can  study  an  isolated  2-DPS.  Hence  the  issue  of  the  vacuum  gap  be¬ 
tween  slabs  periodically  repeated  in  the  surface-normal  direction 
does  not  arise. 

A  single  DMMP  molecule  was  placed  on  the  “top”  surface  of  the 
2-DPS  supercell,  and  a  second  was  placed  symmetrically  on  the 
“bottom”  surface.  This  “two-sided”  approach  preseiwes  the  center 
of  symmetry  of  the  supercell,  thus  avoiding  the  production  of  a  di¬ 
pole  potential,  and  has  been  shown  [43]  to  be  beneficial  in  2-DPS 
studies  of  H2O  adsorption  on  rutile  (110).  The  adsorption  energy 
per  DMMP  is  then  defined  as 

A£ads  =  [£(Ti02  +  nDMMP)  -  £(Ti02)  -  n£(DMMP)  -f  A£(BSSE)]/n 

where  n  =  2,  £(Ti02  +  nDMMP)  is  the  relaxed  total  DFT  energy  with 
adsorbed  DMMP,  £(Ti02)  and  £(DMMP)  are  the  energies  of  bare  Ti02 
and  free  DMMP  and  A£(BSSE)  is  the  counterpoise  correction  for  ba¬ 
sis  set  superposition  error.  Thus  APa^s  <  0  indicates  an  exothermic 
process.  All  technical  parameters  (integration  accuracy,  etc.)  were 
the  same  for  all  quantities  to  be  combined  in  obtaining  A£ads- 
Specifically,  £(DMMP)  was  obtained  in  the  same  manner  as  the 
other  energy  terms;  namely,  by  optimizing  the  geometiy  at  the 
RHF/6-31G(d,p)  level  followed  by  a  single-point  B3LYP/6-31G(d,p) 
calculation.  The  free  DMMP  molecule  was  placed  in  the  lowest- 
energy  gas-phase  conformation  [44,45]  before  starting  the  optimi¬ 
zation  process. 


3.  Results 

3.1.  Relaxation  of  the  bare  surfaces 

Table  1  shows  the  lattice  parameters  obtained  in  the  present 
RHF  calculation  for  bulk  anatase  and,  for  comparison,  bulk  rutile. 
In  both  cases,  the  agreement  with  experiment  [46]  is  seen  to  be 
quite  good.  Table  2  shows  the  relaxation-induced  displacements 
of  atoms  on  and  near  the  bare  unfaceted  surface  in  comparison 
to  previous  RHF  results  [37]  for  a  somewhat  thicker  slab,  together 
with  corresponding  results  [22]  using  DFT  in  the  geometry  relaxa¬ 
tion.  The  agreement  among  the  three  sets  of  results  is  again  quite 


Table  2 

Displacements,  relative  to  the  ideally-terminated  surface,  resulting  from  relaxation  of 
the  unfaceted  anatase  (0  1  0)  surface/ 


Atom'® 

This  work" 

Ref.  [37]'' 

Ref  [22]' 

Layer  1 

Tisc 

0.02;  -0.15 

0 

b 

1 

P 

0.02;  -0.16 

0|| 

0.02;  +0.18 

0.04;  +0.18 

0.04;  +0.18 

Oi 

0.17;  0.00 

0.16;  -0.01 

0.16;  +0.02 

Layer  2 

Tiec 

0.00;  +0.12 

0.00;  +0.11 

0.01;  +0.17 

Oil 

0.00;  +0.07 

0.00;  +0.07 

0.01;  +0.10 

Oi 

0.12;  -0.04 

0.11;  -0.03 

0.13;  -0.03 

^  All  displacements  are  in  Angstroms.  Values  are  given  as  |Ax|;  Az,  where  Ax  is 
the  displacement  in  the  [0  01]  direction,  and  Az  is  the  displacement  along  the 
surface  normal.  A  positive  (negative)  Az  is  outward  (inward).  Layer  1  is  the  surface 
layer.  For  Oj^,  Ax  is  in  a  direction  away  from  the  “channel”  (cf.  Fig.  1).  The  anatase 
(10  0)  and  (0  1  0)  surfaces  are  isostructural. 

0||  refers  to  the  0  atom  bonded  to  three  Ti  atoms  in  the  same  plane  which,  at  the 
surface,  is  an  O3C  site.  0^  refers  to  the  0  atom  with  bonds  to  Ti  in  different  planes 
which,  at  the  surface,  is  an  O2C  site.  See  Fig.  1. 

Unfaceted  (4x1)  2-DPS  model  with  six  layers  after  geometry  relaxation  as 
described  in  the  text. 

^  ( 1 00)-(  1x1)  2-DPS  model  with  eight  layers  relaxed  in  an  RHF  calculation  with 
ECPs  for  Ti  and  0  (see  Ref.  [37]). 

^(10  0)-(l  X  1 )  2-DPS  with  six  Ti  layers  relaxed  in  a  DFT  calculation  using  Car- 
Parrinello  molecular  dynamics  (see  Ref.  [22]). 


Fig.  2.  The  faceted  surface  after  relaxation.  In  correspondence  to  Table  2,  the 
brackets  indicate  the  layer  numbers,  and  the  individual  atoms  are  labeled  by  type. 
For  clarity,  Ti  (0)  atoms  are  labeled  on  the  “top”  (“bottom”)  of  the  slab.  The  Layer  1 
Ti5c(2)  is  obscured  by  an  0  atom.  The  Layer  2  Tisc  site  shown  italicized  in  blue  is  the 
one  used  for  DMMP  adsorption  within  the  facet  (see  text).  The  [0  0  1]  and  [0  1  0] 
vectors  define  the  positive  directions  for  Ax  and  Az,  respectively,  for  the  "top” 
surface  atoms  (cf.  Table  2  and  discussion  in  text).  The  [0  01]  and  [0  1  0]  directions 
are  reversed  on  the  “bottom”  surface. 
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Table  3 

Displacements,  relative  to  the  ideally-terminated  surface,  resulting  from  relaxation  of  the  faceted  anatase  (0  1  0)  surface.^ 


Layer 

Layer  2 

6-Layers 

12-Layers 

6-Layers 

12-Layers 

Atom” 

Ax;  Az 

Ax;  Az 

Atom” 

Ax;  Az 

Ax;  Az 

Tiscd) 

-0.05;  -0.08 

-0.02;  -0.10 

Tiscd) 

+0.06;  -0.25 

0.00;  -0.16 

Tisc  (2) 

-0.01;  -0.13 

+0.02;  -0.10 

Tisc  (2) 

d 

1 

d 

o 

d 

0.00;  -0.16 

0||(1) 

-0.05;  +0.25 

-0.02;  +0.23 

Ti6c(l) 

+0.04;  +0.08 

+0.02;  +0.15 

0||  (2) 

0.00;  +0.21 

+0.03;  +0.23 

Tisc  (2) 

-0.04;  +0.17 

-0.02;  +0.15 

Ox(l) 

+0.13;  0.00 

+0.16;  +0.04 

0||(1) 

+0.05;  +0.09 

+0.01 :  +0.20 

Ox  (2) 

-0.19;  +0.05 

-0.16;  +0.04 

0||(2) 

-0.02;  +0.25 

-0.01 :  +0.20 

0||(3) 

+0.08;  0.00 

+0.04;  +0.06 

0||  (4) 

-0.05;  +0.07 

-0.04;  +0.06 

Oxd) 

+0.24;  +0.01 

+0.19;  -0.02 

Ox  (2) 

-0.19;  0.01 

-0.18;  -0.01 

Ox  (3) 

+0.13;  +0.01 

+0.12;  +0.01 

Ox  (4) 

1 

p 

1 

o 

b 

-0.12;  +0.01 

“  See  Table  2,  footnote  a.  “6  Layers"  refers  to  the  2-DPS  with  a  (4  x  2)  supercell  and  6- 
supercell  and  12-layers  (before  faceting).  See  Section  3.1  for  further  details. 

See  Table  2,  footnote  b,  and  Fig.  2. 

good.  The  relaxation  is  essentially  identical  for  all  equivalent  atoms 
in  the  present  (4x1)  cell,  which  indicates  no  tendency  toward 
reconstruction.  Fig.  2  and  Table  3  show  corresponding  results  for 
the  faceted  surface,  which  are  more  complex.  For  the  2-DPS  de¬ 
scribed  above,  with  6-layers  before  faceting,  there  are  differences 
in  the  displacements  of  equivalent  atoms  having  different  posi¬ 
tions  in  the  [001]  direction.  For  example,  the  two  Layer  1  Tisc  sites 
should  be  equivalent  (cf.  Fig.  lb)  but  exhibit  somewhat  different 
displacements  during  relaxation.  However,  equivalent  atoms  at 
different  positions  in  the  [10  0]  direction  have  essentially  identical 
relaxations. 

These  anomalies  result  from  the  fact  that  a  thickness  of  only  6- 
layers  is  insufficient  for  complete  relaxation  when  the  slab  is 
thinned  by  faceting.  This  is  shown  by  results  for  a  slab  with  12- 
layers  before  faceting  and  a  (1  x  2)  supercell.  Reducing  the  cell 
dimension  in  the  [1  0  0]  direction  is  permissible  in  view  of  the  re¬ 
sults  obtained  for  the  6-layer  slab.  For  the  thicker  slab  the  dis¬ 
placements  of  equivalent  atoms  are  now  essentially  identical, 
but  these  are,  in  most  cases,  not  very  different  from  those  for 
the  thinner  slab.  The  signs  given  for  Ax  in  Table  3  are  significant. 
Thus  Ox(l)  and  Ox(2)  in  the  Layer  1  displace  away  from  the  facet, 
and  Ox(l)  and  Ox(2)  in  the  Layer  2  displace  away  from  the  “chan¬ 
nel”.  In  Layer  6  of  the  12-layer  slab  (not  shown)  the  displacements 
from  the  bulk  positions  are  all  <0.03  A  in  magnitude,  indicating 
reasonably  good  convergence.  Comparison  of  Figs.  1  and  2  shows 
that  relaxation  leads  to  an  atomic-scale  roughening  of  the  surface, 
as  noted  previously  [22,37],  with  0||(Ti5c)  displacing  outward  (in¬ 
ward)  along  the  surface  normal.  Finally,  the  density  of  states  for 
either  surface  (not  shown)  revealed  no  surfaces  states  in  the  band 
gap,  in  agreement  with  previous  results  [10,17]  for  the  relaxed 
unfaceted  surface. 

3.2.  Molecular  adsorption  of  DMMP 

3.2.1.  The  unfaceted  surface 

Fig.  3  shows  the  relaxed  configuration  of  DMMP  molecularly 
adsorbed  on  either  surface.  For  the  unfaceted  surface  (Fig.  3a) 
AEads  = -21.0  kcal/mol  was  found  after  a  BSSE  correction  of 
+6.5  kcal/mol.  The  Ti-  ■  OP  and  TiO=P  bond  lengths  are  2.09  and 
1.50  A  respectively,  and  the  Ti  -  0=P  bond  angle  is  138°.  There 
appear  to  be  two  C— H-  ■  ■02c  interactions  that  meet  the  H— bond 
criteria  [47]  for  interatomic  distances  and  angles.  For  Bond  1  (2) 
the  H  -  0  and  C  -  0  distances  are  2.28  and  3.29  A  (2.65  and 
3.67  A)  respectively,  and  the  C— H-  ■  O  angle  in  either  case  is 
157°.  Other  possible  C— H  -  0  interactions  (based  on  the  H  -  0 
distances)  involve  much  smaller  angles  generally  indicative  of  a 


-layers  (before  introduction  of  the  facet).  "12  Layers”  refers  to  the  2-DPS  with  a  (1  x  2) 


weaker  bond.  It  is  noted  that  AEads  for  the  (0  1  0)  surface  is  close 
to  2-DPS  results  obtained  elsewhere  [42]  for  molecular  adsorption 
of  DMMP  on  the  anatase  (1  0  1)  and  rutile  (1  1  0)  surfaces  (-19.1 
and  -17.6  kcal/mol,  respectively).  The  small  variation  in  A£ads  re¬ 
flects  the  trend  in  surface  energies  given  above,  with  AEads 
decreasing  in  magnitude  with  increasing  surface  stability.  On  all 
surfaces,  adsorption  reverses  the  inward  displacement  of  the  Tisc 
bonding  site  that  occurs  when  the  bare,  ideally-terminated  surface 
relaxes.  This  can  be  seen  clearly  in  Fig.  3a  as  an  outward  displace¬ 
ment  of  the  Tisc  adsorption  site  relative  to  the  neighboring  Tisc 
sites. 


o 


Fig.  3.  Sections  of  the  (a)  unfaceted  and  (b)  faceted  surface  after  relaxation  with 
adsorbed  DMMP  in  the  lowest-energy  conformations.  For  clarity,  only  sections  of 
the  entire  slab  supercells  are  shown.  The  dashed  lines  show  H-bonding  interactions 
identified  as  such  on  the  basis  of  the  02c-  ■  -H  and  02c-  ■  -C  distances  and  the 
02c-  ■  -H— C  angle  (see  text).  In  (b)  the  adsorption  site  is  the  one  designated  Layer  2 
Ti5c(l )  in  Fig.  2.  Note  that  the  orientation  in  (b)  is  reversed  relative  to  that  in  Fig.  2. 
The  green,  black  and  gray  spheres  represent  P,  C  and  H  respectively.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is 
referred  to  the  web  version  of  this  article.) 
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Other  configurations,  differing  from  that  in  Fig.  3a  by  rotation  of 
the  DMMP  by  90°  or  180°  about  the  P=0  bond,  were  also  consid¬ 
ered  as  starting  points  in  the  geometry  optimization.  These  differ 
as  to  which  combinations  of  CH3  groups  are  positioned  favorably 
for  C— H-  ■  ■O2C  bond  formation.  At  the  RHF  level  used  in  geometry 
optimization  the  90°  (180°)  rotations  relaxed  to  energies  which 
were  higher  by  about  2.3  (9.0)  kcal/mol  than  that  of  Fig.  3a.  An¬ 
other  starting  geometry,  differing  from  that  in  Fig.  3a  by  a  rotation 
of  the  DMMP  by  180°  about  the  Ti-  ■  OP  bond,  was  also  considered. 
This  structure  makes  use  of  the  slightly  more  distant  O2C  sites  for 
possible  H-bonding.  Flowever,  the  relaxed  structure  was  less  stable 
by  ~15  kcal/mol  at  the  RHF  level.  Still  other  rotations,  by  90°  about 
the  Ti-  ■  OP  bond  or  by  180°  about  both  the  P=0  and  Ti-  -  -OP  bonds 
also  relaxed  to  structures  which  were  unstable  relative  to  that  in 
Fig.  3a. 

The  fact  that  none  of  these  starting  structures  relaxed  to  the 
apparently  most  stable  configuration  (Fig.  3a),  but  instead  re¬ 
mained  in  approximately  the  initial  conformation,  indicates  that 
there  is  a  significant  barrier  to  the  bond  rotations  needed  to 
accomplish  the  transformation.  This  may  arise  from  the  need  to 
break  one  or  more  C— H-  -  -O2C  bonds  in  order  to  accomplish  the 
rearrangement.  The  existence  of  this  barrier  made  it  necessary  to 
consider  several  different  initial  structures  in  the  geometry 
optimization. 

Finally,  dative  bonding  between  a  Tisc  and  a  methoxy  O  atom  is 
also  a  possible  mode  of  adsorption.  This  has  been  considered  else¬ 
where  [42]  for  the  rutile  (1  1  0)  surface  and  been  found  to  give  a 
negligible  AFads-  In  the  present  case,  an  attempt  was  made  to  ad¬ 
sorb  DMMP  on  the  unfaceted  surface  via  a  Tis^-  -  -O  bond  to  a  meth¬ 
oxy  O  atom.  At  the  RHF  level  used  in  geometry  optimization  the 
resulting  structure  was  less  stable  than  that  in  Fig.  3a  by 
~25  kcal/mol,  indicating  an  unstable  configuration.  Likewise,  an 
initial  structure  with  a  Tis^-  -  -0=P  bond  together  with  a  Tis^-  -  -0 
bond  to  a  methoxy  at  an  adjacent  Tisc  site  relaxed  to  one  with  only 
the  Tisc-  -  -0=P  bond. 

3.2.2.  The  faceted  surface 

For  the  faceted  surface  (Fig.  3b),  DMMP  was  placed  symmetri¬ 
cally  on  both  slab  faces  at  the  Ti5c(l)  site  of  Layer  2  shown  in 
Fig.  2.  The  slab  at  this  point  is  four  layers  thick,  which  should  suf¬ 
fice  for  a  nearly-converged  Afads.  and  is  preferable  to  the  Layer  2 
Ti5c(2)  site  where  the  slab  is  five  layers  thick.  Previous  work  on 
the  adsorption  of  CO  [48]  or  H2O  [49,50]  on  rutile  (1  1  0)  shows 
that  a  thin  2-DPS  with  an  even  (odd)  number  of  Ti  layers  tends 
to  slightly  underestimate  (significantly  overestimate)  the  fully- 
converged  magnitude  of  AFads-  Similar  odd/even  oscillations  in 
the  (1  0  0)  surface  energy  are  seen  for  thin  anatase  2-DPS  slabs 
[37[. 

For  the  faceted  surface,  AFads  =  -21.3  kcal/mol  was  found  using 
the  same  AF(BSSE)  as  above.  The  Ti-  -  -OP  and  TiO=P  bond  lengths 
are  2.04  and  1.49  A  respectively,  and  the  Ti-  -  -0=P  angle  is  140°.  All 
results  are  close  to  those  given  above  for  adsorption  on  the  unfac¬ 
eted  surface.  The  dashed  line  indicates  one  of  several  possible  H- 
bonding  interactions,  based  on  the  H-  -  -0  distances  which  are  in 
the  range  of  2.39-3.00  A.  The  indicated  bond  is  believed  to  be 
the  strongest,  with  0-  -  -H  and  C-  -  -0  distances  of  2.39  and  3.29  A 
respectively  and  a  C— H-  -  -O  angle  of  140°.  Other  interactions  have 
smaller  C— H-  -  -0  angles  and/or  longer  H-  -  -O  distances  and  are 
therefore  probably  weaker  [47].  It  is  possible,  however,  that  the 
cumulative  effect  of  several  weak  H-bonds  might  make  a  signifi¬ 
cant  contribution  to  the  total  AFads.  A  quantitative  discussion  of 
H-bonding  requires  the  characterization  of  bond  critical  points 
using  the  atoms  in  molecules  theory  [51]  which  is  not  available 
at  present. 

Alternative  starting  structures,  differing  from  that  in  Fig.  3b  by 
rotations  of  70°  or  140°  about  the  P=0  bond,  were  also  considered 


in  order  to  position  various  CH3  groups  for  possible  H-bonding  to 
O2C  sites  in  Layers  1  and  2.  However,  these  were  less  stable  by 
about  3.5  and  7.2  kcal/mol  respectively  at  the  RHF  level  used  in 
geometry  optimization.  As  in  the  case  of  the  unfaceted  surface,  nei¬ 
ther  of  these  initial  structures  relaxed  to  the  most  stable  configura¬ 
tion  during  geometry  optimization  which  suggests  a  significant 
barrier  to  the  necessaiy  bond  rotation.  Large  rotations  about  the 
Ti-  -  -OP  bond  were  also  considered  but  were  hindered  by  steric 
interactions  between  the  CH3O  groups  and  the  “step  edge”  be¬ 
tween  Layers  1  and  2. 

As  noted  above,  both  Tisc  sites  in  Layer  2  (Fig.  2)  are  equivalent 
on  the  semi-infinite  surface  but  can  differ  for  a  thin  2-DPS.  Adsorp¬ 
tion  at  the  Ti5c(2)  site  of  Layer  2  was  evaluated  for  DMMP  in  the 
same  configuration  as  at  the  Ti5c(l)  shown  in  Fig.  3b.  The  bond 
lengths  and  angles  in  the  relaxed  structure,  including  those  of 
the  putative  dominant  H-bond,  were  virtually  identical  to  those 
found  for  adsorption  at  Ti5c(l),  but  AF^ds  was  larger  in  magnitude 
by  ~4.3  kcal/mol  at  the  RHF  level  used  in  geometry  optimization. 
On  the  basis  of  the  above  discussion,  this  is  assumed  to  represent 
an  overestimate.  Adsorption  outside  the  facet  (i.e.,  at  a  Layer  1 
Tisc  site  in  Fig.  3b)  was  compared  with  adsorption  inside  the  facet 
(i.e.,  at  the  Layer  2  Tisc  site  as  shown  in  Fig.  3b).  The  starting  geom¬ 
etry  was  the  same  as  those  identified  above  as  the  most  favorable 
and  shown  in  Fig.  3a  and  b.  At  the  RHF  level  the  resulting  AFads  was 
3.9  kcal/mol  smaller  in  magnitude  than  for  adsorption  inside  the 
facet.  The  overall  conclusion  is  that  the  results  do  not  depend 
greatly  on  the  choice  of  adsorption  site  but  that  the  structure  in 
Fig.  3  b  is  the  most  favorable  for  molecular  adsorption  on  the  fac¬ 
eted  surface. 

3.3.  Dissociative  adsorption  of  DMMP 

Fig.  4  shows  models  for  a  possible  first  step  in  the  dissociation 
of  DMMP  on  surfaces  with  and  without  faceting.  The  structure  is 
obtained  by  breaking  a  PO— CH3  bond,  moving  the  CH3  to  an  O2C 
site  and  then  forming  a  (CH3)(CH30)P(— O— )2  bridge  between 
two  Tisc  sites.  This  is  only  a  formal  description  and  is  not  intended 
to  represent  the  actual  reaction  mechanism.  The  end  product  is 
consistent  with  the  results  of  X-ray  photoemission  spectroscopy 
(XPS)  and  temperature-programmed  desorption  (TPD)  experi¬ 
ments  for  DMMP  adsorbed  at  room  temperature  on  OH-free  sin¬ 
gle-crystal  rutile  (110)  surfaces  prepared  and  maintained  in 
UHV  [29,30].  These  data  show  mainly  molecular  adsorption.  The 
relatively  small  amount  of  decomposition  that  does  occur  leads 
to  the  desorption  of  CH3,  but  not  CH3O,  in  TPD.  The  yield  of  CH3 
is  distinct  from  the  cracking  of  desorbed  DMMP  itself  and  is  consis¬ 
tent  with  CH3  adsorbed  at  a  bridging  O2C  site  as  opposed  to  a  ter¬ 
minal  CH3— O— Ti  site.  Desorption  of  CH3O  from  such  a  terminal 
site  would  reform  a  Tisci  whereas,  desorption  of  CH3  would  yield 
an  unstable,  doubly-unsaturated  0  atom  unless  a  Tisc  site  were 
available  for  forming  an  O2C  bridge.  However,  these  data  were  ob¬ 
tained  for  rutile  (1  1  0)  and  might  not  necessarily  apply  directly  to 
anatase  (0  1  0). 

A  thorough  examination  of  all  a  priori  possible  dissociation 
products  was  not  performed  since  the  goal  here  is  to  determine 
the  extent  to  which  faceting  affects  dissociation,  and  the  relative 
energies  for  molecular  vs.  dissociative  adsorption  for  a  plausible 
structure.  Several  other  dissociation  structures  are  considered  else¬ 
where  [42]  for  DMMP  on  rutile  (1  1  0)  and  anatase  (101).  The 
starting  configuration  was  based  on  the  most-favorable  structures 
for  molecular  adsorption  discussed  above  and  shown  in  Fig.  3.  For 
the  unfaceted  surface,  one  alternative  structure  was  considered  in 
which  the  CH3  was  placed  on  one  of  the  two  equivalent  O2C  sites 
which  are  closest  to  the  bridge  (i.e.,  one  step  to  the  right  of  that 
occupied  in  Fig.  4b).  This  alternative  site  is  marked  by  an  arrow 
in  Fig.  4b.  The  OCH3  on  the  bridge  was  also  rotated  90°  about  the 
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Fig.  4.  Similar  to  Fig.  3  but  showing  structures  considered  for  the  first  step  in  DMMP  dissociation.  The  upper  panel  shows  the  unfaceted  surface  viewed  along  the  (a)  [1  0  0| 
and  (b)  [0  0  1]  direction,  (c)  Shows  the  faceted  surface  viewed  along  the  [1  0  0]  direction.  Possible  H-bonding  interactions  have  been  omitted  for  clarity.  The  arrow  in  (b) 
marks  an  alternative  location  for  the  CH3  fragment  which  is  discussed  in  Section  3.3. 


P— OCH3  bond  so  as  to  move  the  CH3  away  from  the  O2C— CH3,  thus 
minimizing  steric  repulsion  between  the  two  CH3  groups  and 
allowing  C— H-  ■  O  bonding  between  the  O2C— CH3  and  the  0  atom 
in  the  P— O— CH3  group.  At  the  RHF  level  the  geometiy-optimized 
A£ads  was  negligibly  different  from  (~1  kcal/mol  larger  than)  that 
for  the  structure  in  Fig.  4b. 

The  BSSE-corrected  AEaUs  is  -32.7  and  -22.7  kcal/mol,  respec¬ 
tively,  for  the  unfaceted  and  for  the  faceted  surface.  These  are  lar¬ 
ger  than  the  energies  given  above  for  molecular  adsorption  but 
only  slightly  so  for  the  faceted  surface.  This  indicates  that  dissoci¬ 
ation  is  thermodynamically  favored  to  a  degree  which  depends  on 
faceting.  As  noted  above,  faceting  exposes  the  more  stable  (Oil) 
surface,  which  is  isostructural  with  the  (1  01).  Consistent  with 
this,  a  2-DPS  calculation  done  elsewhere  [42]  for  the  dissociative 
adsorption  of  DMMP  on  the  anatase  (10  1)  surface,  in  a  structure 
like  that  studied  here,  gave  a  BSSE-corrected  AEads  of  -25.0  kcal/ 
mol  which  is  close  to  the  present  result  for  the  faceted  (0  1  0). 
Although  no  relevant  data  are  available  for  completely  OH-free 
anatase,  DMMP  adsorption  on  OH-free  rutile  (1  1  0)  at  room  tem¬ 
perature  is  largely  molecular  [29,30],  If  the  same  is  true  for  anatase 
(0  1  0)  then  the  present  results  suggest  that  dissociation  under 
these  conditions  is  kinetically  limited.  This  is  consistent  with 
experimental  results  for  partially-hydroxylated  Ti02  powder 
[31,32]  which  indicate  that  surface  OH  is  directly  involved  in  the 
decomposition  of  molecular  DMMP,  possibly  via  a  nucleophilic  at¬ 
tack  on  the  P  atom  by  the  0  atom  of  a  Ti— OH  group,  and  that  this 
can  occur  below  room  temperature.  In  the  present  interpretation 
the  Ti— OH  is  necessaiy  to  overcome  the  kinetic  barrier  to 
decomposition. 

4.  Conclusions 

This  is,  to  our  knowledge,  the  first  computational  study  of  the 
faceted  anatase  (0  1  0)  (or,  equivalently,  (1  0  0))  surface  either  with 
or  without  an  adsorbate.  On  the  bare  surfaces,  similar  relaxations 
occur  for  corresponding  atoms  with  and  without  faceting.  These 
consist  mainly  of  inward  (outward)  displacements  of  Tisc  (Tisc) 
atoms  and  outward  displacements  of  O3C  atoms.  Surface  O2C  atoms 


displace  mainly  parallel  to  the  surface.  The  energies  and  structures 
for  molecular  adsorption  of  DMMP  are  similar  in  either  case  and 
are  also  similar  to  those  [42]  for  rutile  (110)  and  anatase  (101) 
surfaces.  Adsorption  occurs  preferentially  via  a  Tisc-  ■  0=?  dative 
bond  with  evidence  of  C— H-  ■  0  bonding  between  methoxy  groups 
and  O2C  atoms.  The  results  indicate  that  molecular  adsorption  of 
DMMP  on  HSA  Ti02  powders,  for  which  the  anatase  (0  1  0)  surface 
is  a  major  constituent,  can  be  understood  without  necessarily  tak¬ 
ing  account  explicitly  of  the  faceting.  As  a  corollary,  the  results  also 
suggest  that  data  for  molecular  adsorption  of  DMMP  on  the  easily- 
obtained  rutile  (110)  surface  should  be  transferable  to  the  less 
common  anatase  (0  1  0)  surface.  However,  this  applies  only  to 
the  pristine  surfaces  considered  here  which  are  free  of  defects 
and  adsorbed  OH.  The  situation  may  change  when  such  imperfec¬ 
tions  are  introduced. 

A  possible  first  step  in  the  dissociative  adsorption  of  DMMP  on 
OH-free  surfaces  was  defined  on  the  basis  of  experimental  data 
[29,30]  for  OH-free  rutile  (1  1  0)  single-crystal  surfaces.  On  the 
non-defective  and  OH-free  surface  the  product  consists  of  an 
O2C— CH3  bond  at  a  Ti— O2C— Ti  site  and  a  (CH3)(CH30)P(— 0— [2 
bridge  between  two  Tisc  sites.  AEa^s  for  this  process  indicates  that 
it  is  thermodynamically  favored  over  molecular  adsorption  but  to  a 
degree  that  depends  on  faceting.  Faceting  exposes  the  more  stable 
(011)  (or,  equivalently,  (101))  surface  and  reduces  the  magni¬ 
tude  of  AEads  for  dissociative  adsorption,  shifting  it  toward  the  re¬ 
sult  obtained  for  the  actual  (101)  surface.  Experiments  on  fully 
OH-free  anatase  (0  1  0)  surfaces  are  needed  in  order  to  determine 
whether  dissociation  is  kinetically  accessible  without  the  interven¬ 
tion  of  Ti— OH  sites. 

Note  added  in  proof 

A  few  important  references  were  overlooked  in  the  original 
manuscript.  Barnard  et  al.  [52-54]  have  performed  ab  initio  DFT 
studies  of  the  anatase  (1  0  0),  among  several  other  Ti02  surfaces. 
The  structures  of  the  clean  surfaces  and  the  effects  of  adsorbed  H 
atoms,  0  atoms  and  H2O  were  analyzed.  Ignatchenko  et  al.  [55] 
also  used  ab  initio  DFT  to  study  the  molecular  and  dissociative 
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adsorption  of  H2O  on  anatase  (10  0)  and  other  Ti02  surfaces.  Gong 
and  Selloni  [56]  performed  DFT  studies  of  the  adsorption  of  H2O, 
CH3OH  and  HCOOH  on  the  anatase  (10  0)  surface  and  on  steps 
on  the  (10  1)  with  a  structure  similar  to  that  of  the  extended 
(1  0  0)  surface. 
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